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Abstract 

We report a very detailed relation about the study, by a con- 
strained Monte Carlo dynamic of a 4D EA spin glass model (J = ±1). 
In particular we concentrate our attention on the study of the be- 
haviour of the system under different dynamical parameters in order 
to optimize our analysis and try to understand the low energy states 
structure. We find that in the thermodynamic limit this structure 
assumes just ultrametric features as already established for the SK 
model in the SRSB (Spontaneous Replica Simmetry Breaking) the- 
ory. 

Introduction 

At present there are two different and contrasting approaches about the de- 
scription of the spin-glass phase of a finite-dimensional EA model. The first, 
the rephca approach, follows the solutions of the SK model describing the 
nature of spin-glass phase at finite dimensions similar to the one of the mean- 
field theory The second approach is based on the phenomenological 
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theory of droplets and predicts a spin-glass phase dominated by only one 
equilibrium state [0. The difficulties in the research of analytic results in 
this field imposes that a large part of the work is delegated to the numerical 
simulations. A large number of works (||19[|10[ 131 0) supports the replica 
approach which at present seems the more appropriate for the description of 
realistic spin-glasses. 

In this work we study the four-dimensional Ising Spin Glass at zero exter- 
nal magnetic field by means of Monte Carlo numerical simulations for several 
small sizes N = with L=3, 4, 5, 6, 7, 8 with coupling quenced ±J. 

In detail we will concentrate our attention on the structure of low energy 
states of the model by a constrained dynamic that will be here examinated 
very in detail (see also [|T|). In the ffist and second section we will briefly 
mention the analytic results obtained with the SK model and we will add 
some comments about the main numerical works undertaken about the ultra- 
metric problem in spin-glasses (see also 0). In section 3, we will introduce 
the metric used for this model. Sections 4, 5, and 6 are dedicated to the 
description of the simulated system and to the Monte Carlo constrained dy- 
namic realized. Finally the last sections report the results of the numerous 
simulations effected. 



1 The ultrametricity in the mean field theory 

Rigorous analytic studies (1T^0|§[) have clearly demonstrated the ultramet- 
ric structure of low energy states of the SK model. The starting point of this 
analysis is the calculation of the probability for three pure states ai, 02,(^3 
to have mutual overlaps 

gi = g"2"3^ g2 = g"^"S g3 = g"^"^ . (1) 

One can demonstrate that subsists the relation 

PiQu 92, gs) = lP{qi)x{qi)6{qi - q2)S{qi - ^3) + 

+ ^{P{qi)P{q2)0{qi - q2)S{q2 - ga) + 2 permutations } . (2) 

The study of this equation |^ tells us that P{qi,q2,q3) is null except when 
two of the overlaps are equals and not larger of the third one. If we define 
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the distance between two pure states as 

dl, = ]^T.i<-^^)\ (3) 

i 

(this is not the only way), can be easily verified that 

dip = ^Qea -qafs) ■ (4) 

In this contest the equation (H) establishes that triangles builded with three 
configurations, chosen in agreement with their Boltzmann weights, are always 
either equilaters or isosceles, and in this last case the different side is the 
smallest. In this prospective, we need to replace the triangular inequality 

dab < dac + dbc (5) 

with the more restrictive condition 

dab < max{dac,dbc) ■ (6) 
A space with this metric is called ultrametric. 



2 Ultramericity in literature 

Attempts of analysis with calculator of ultrametric properties of finite di- 
mensional spin glasses are very few in literature ( ) and have not given 
any decisive contribute. The more remarkable results of interest for our work 
are probably those obtained by Ciria, Parisi and Ritort in 4 dimensions 
at T=1.4. Now we discuss the crucial elements of the philosophy followed 
in that work. For every triad of three overlaps that one can build (see also 
section 4) from three different replicas of the system (with the same disor- 
der realization but an autonomous dynamic) one analyses, for a fixed bigger 
overlap in an established range of values, the difference between the middle 
overlap and the smallest excluding terms where the triangular inequality is 
violated from finite size corrections. For this purpose one defines for every 
overlap a distance 

dal3 = [2{qma. - q""^)]'^' , (7) 
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where Qmax is the largest overlap relative to the thermodynamic limit 

Qmax = lim qi^^ , (8) 

L—>oo 

(which can be obtained by the scaling relation P{q > qmax) ~ -Pmaa;^(A/'(n — 
p 

■'-'■$H§) ^^'i '^'^^ consider only the overlap triads that satisfy the in- 

equality 

d'max ^ djYiin ~l~ dYniddle • (9) 

The results of this analysis have not given a decisive answer to the problem. 
In effect if we observe figure (0), where is indicated in abscissas the value of 
^1 = Qmedio — Qmin and in ordinates P{Sq), we note that a dependence on the 
curves from the volume is not determined with sufficient statistic precision. 
By these simulation people thought that it were possible to increase the 
volume size to obtain a precise indication of a ultrametric behaviour in the 



infinite volume limit. Similar works about the SK model |T6[|T^ have given 
very satisfying results pointing out clearly the ultrametric structure of the 
states. The first of these works is from Bhatt Young |T6[. The result is shown 



in figure @. In this case the curve tends explicitly to narrow increasing the 
lattice volume showing the existence of the predicted structure between the 
states. 



3 The distance 

There are two big problems for a direct analysis of the ultrametric structure: 

1) There are corrections for small N that are not known analytically 

2) The triangular inequality can impose some bonds that are not easily 
distinguishable from the more restrictive ones provided from the ultrametric 
inequality. 

To understand the last affirmation we take a triangle with sides di, d2, 
d^ and let d^ the smallest one. The triangular inequality imposes that 

<rfi <c?2 + 4 (10) 

to compare with the ultrametric request 

d^ = d2. (11) 
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A correction due to the limited value of N causes some violations in the last 
equation 

di = d2± 0(l/iV^) . (12) 

If (is is fairly little, the imposed conditions from the two inequality are prac- 
tically indistinguishable. To get over these problems it is necessary to define 
a metric which makes the two conditions as distinct as possible. 

We define (see also ^) the square distance between the two pure states 
with overlap q as 

,2 QEA-q 

d = — , (13) 

so when q = —qEA "we will have c?^ = 1, the lowest overlap imposes a largest 
square distance, while for q = qEA we have that d'^ = 0, the largest overlap 
involves a lowest square distance. The triangular inequality imposes that 
when we take three pure states the mutual distances satisfy the relation 

dl3 + ^23 > di2 . (14) 

Squaring both members we have 

idi3 + d^^f > dl^ (15) 

and so 



{qEA-qiz) ^ (ggA-923) ^ ry VOEA- qisVqEA- q23 ^ qEA- qi2 ^^g^ 

2g£;A "iqEA '2qEA ~ '^qEA 

We fix now two of the three overlaps, for exemple 

gi3 = 923 = g^^^ (17) 

(we take for convenience of notation di2 = d and qi2 = q)- 
From the equation (^) we have 

^{qEA - q^n >qEA-q>0, (18) 

because when di^ > ^23 then the triangular inequality imposes 

di2 > (rfi3 - ^23) ■ (19) 
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In this way the overlap q (the observable we will study in our Monte Carlo 
simulations) is forced from the triangular inequality to satisfy the relation 



qEA>q> 4g^'" - 'iqEA • (20) 
It is reasonable to choose for our simulations for exemple 

g^'" = \qEA ■ (21) 
5 

With these conditions the triangular inequality imposes that 

- ^-qEA <q<qEA (DT bound) , (22) 
while the ultrametric one imposes 

\qEA<q< qEA (UM bound) . (23) 
5 

Now the difference between these two situations is evident. When N is 
limited (and little) these inequalities will be violated. The possible ultra- 
metric structure will completely manifest itself only in the limit of infinite 
volume. We will try to understand if the more we increase the lattice size 
the better the ultrametric bound is satisfied. 

We will use for qea the thermodynamic limit calculated in a work of 
Ciria, Parisi and Ritort |12|; at T = 1.4, qea = -54. So we will take for our 
simulations 

g^*^ ~ .21 (24) 



4 The model 

The studied model in this work is a Ising Spin Glass in four dimensions in 
absence of external magnetic field. It is defined within a 4D lattice, where the 
lattice sites are individualized from the integer values of the fourdimensional 
vector i. In every site is defined a Ising spin (cTj = ±1). The couplings Jjj 
are distributed in binary way 

nJ^,) = l^iJ,, - 1) + ^5iJ,, + 1) (25) 



6 



so that they assume the vahies Jij = ±1 with the same probabiUty and their 
action range is hmited to only first neighbouring spins. The hamiltonian of 
the system is 

H[a] = -WJ^,^^^^j ■ (26) 



2.. 



Prom these definitions, if we consider the overlap between two thermalized 
replicas of the system 



1 ^ 



Q^^Ea'^^ ' (27) 

i=i 

it is possible to calculate easily the overlaps distribution function 



1 ^ 



Piq) = PAq) = {S{q - ]v Ep^^O^ (28) 



which is the observable upon which we will concentrate our attention. 



5 The dynamic 

The exploration of the phase space at a particular temperature {T—IA in 
our case) is directed principally by two algorhythms 

• Monte Carlo Algorhythm 

• Bound Algorhythm 

The Monte Carlo Algorhythm used in our simulations generates a single- 
spin-fiip Glauber dynamic, while the Bound Algorhythm, the very heart of 
the simulations, acts re-examining the Monte Carlo flips on the ground of 
a bound imposed on the overlaps. They work in this way: One takes three 
replicas of the system, 5*1, 5*2, 6*3, with the same realization of the disorder 
but different initial spin configurations. One builds the overlaps 

N ^ N ^ N 

1m = J^J1^1^2 Q{1S) = J^Jl^l^S 9(23) = ^ E '^2'S'3 (29) 

and fixes the value of gf(i3) and 5(23) equal to q^'^^ (as before estabhsched) and 
defines the tolerance parameter inside of which one consents same fluctua- 
tions 

9(13) = 9(23) = [9^*" ± e] . (30) 
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The spin changes proposed from the Monte Carlo Algorhythm that satisfy 
this request will be accepted. The others will be rejected. 
Now see in detail how to operate the chosen. Suppose that in a particular 
instant at a particular site i the Monte Carlo Algorhythm proposes to the 
Bound Algorhythm the spin changes 

Sl~^Sl S\^Sl Sl^Sl. (31) 

Let be 5(^3) = j^T.^,S\Sl and 5^3) = j^T^tiS^Sl the overlaps 
calculated from the previous accepted configurations. Consider to begin the 
spin change S\ ^ SI. One calculates the product = SlSl and tries to 
replace = S{Sl with 5^3' in ^(^3^ 

QiTs) = qU) - lis + Q\3- (32) 
Only in the case in which qj0j e [q^''^ — s, q^^^ + e] we accept the configuration 

Consider now the spin change 5*2 Si and one imposes an argument similar 
to the previous for the overlap q^2i)- ^^^^ ™ ^^^^ accept the change 

only if g^23) ^ ~ Q.^^^ + ^l- 

One examines to finish the spin change S^. We calculate q\^' = 

S^Sl e ^23" = 'S'3'S'2 and replace ^5^3 with ql^" in 5^3^ and q^^ with 

i " ■ V 
?23 m ?(23) 

?(i3) = qIs) - Qls + Qls" (33) 

9(23) = - 923 + Q23' ■ (34) 

Only if both of them will be included in the defined interval we will accept 
the replacement SI ^ S^. 

For our simulations we have estabhshed after several tests that the better 
value of the tolerance parameter is £ = 0.04 (see section 9). 

6 Initial configurations 

The initial spins configurations are so arranged (it is not the only possible 
way) : 
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1) One generates in a random way the spin configuration of tlie replica 5*1. 

2) One builds the replica S3 so that 60% of spins are (chosen in a random 
way) equals to the correspondent spins of the replica Si and the others are 
opposites. 

3) The replica ^2 will be builded in a way such that 40% of spins (chosen 
casually) are opposites to those correspondent of the replica S3 and the others 
are equals. 

These choices allow to begin the simulations in a random spin situation and 
from values of constrained overlaps 

9(23) = g(i3) - 0.2 . (35) 
g(i2) varies with the statistic shows in figure (|^). 



7 Free and constrained P(q) 

The first part of this work regards the study of the overlaps distribution P(q) 
at the temperature T=1.4. This work is useful for two aims. The first is to 
find some good annealing schedules for all lattice sizes (taken successively as 
foundation for the study of the ultrametricity) while the second regards the 
understanding of the behaviour of this distribution for g ~ 0, that is to say 
P(0). The analysis we have done confirms the results already obtained in 
other works [0[|14|; it seems that there is a clean indipendence of P(0) from 
the size of system in agree with the Parisi's theory. This situation seems to 
confirm the existence of a large number of pure states in the glass-phase that 
are describable with the theory of SRSB (Spontaneous Replica Symmetry 
Breaking). In the figure (^) we report the results of these simulations. For 
lattices of size L=3 and L=4 the averages over the disorder are effected with 
1200 different realizations, for L=5 with 700 samples and for the lattice of 
size L=6 we have averaged with 350 samples. 

As already described in detail, the study of the ultrametricity are made 
with a constrained dynamic. We have fixed two of three overlaps built with 
three replicas of the system and we have observed the behaviour of the not 
constrained one changing the lattice size. The triangular inequality , consid- 
ering the bounds imposed to the dynamic, tell us that 

7 

- -Qea <q< Qea (DT bound) (36) 
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while the ultrametric one 



-qEA <q<qEA (UM bound) (37) 

(?EA=0.54). For finite (and httle) N these inequahties will be violated. A 
possible ultrametric structure will reveal oneself completely only in the ther- 
modynamic limit. The annealing schedule used for these simulations contains 
larger numbers of Monte Carlo Steps compared with the ones used for the 
free P(q) thermalization. Not only the number of MCS is on average larger 
that the previous case, but for sizes L=5, L=6, L=7 and L=8, the step of 
the simulated annealing is reduced to one half. Of course the last operation 
caused a greater precision in the research of the equilibrium, but so the sim- 
ulation times are resulted clearly longer. This choice (extremely prudent) is 
operated because, cause the originality of the dynamic, it is not still defined 
any thermalization method (however the annealing schedule relative to the 
thermalization of the free P(q) would be more than enough to assure the 
equilibrium of our system). 

Figure (|^) shows the result of simulations varying of the lattice size. 

The behaviour of this curves point out a clear dependence of constrained 
P(q) from the volume of system. Increasing the lattice side we observe that 
the function tends to assume such a shape that the ultrametric bound is more 
and more satisfied. Figure (^ shows, with vertical lines, both triangular 
bound and 

- 0.75 <q< 0.54 (DT bound) , (38) 

ultrametric one 

0.21 <q< 0.54 (UM bound). (39) 

It can be observed that increasing lattice volume there is a systematic 
shift of the peak toward overlap values consented from the ultrametric bound. 
Similar remarks can be made about tails relative to negative overlaps because 
they too tend toward values indicated from vertical continuous lines. So it is 
possible to deduce a clear tendency for large volumes to a fully ultrametric 
behaviour of the system. The study of the variation of the largest overlap 
ilmax) "with the lattice size [Q (see figure consents moreover to establish a 
limit to the peak shift of constrained P(q). Best fits made with the function 
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=n°° +— (40) 

indicate that g^^^ = 0.31 ±0.09. It is to say that in the thermodynamic hmit 
the constrained P(q) will assume her largest value just inside the ultrametric 
bound. 

Further confirmation of the ultrametric behaviour of the system can be 
seen studying the variation of integral (see also []T|) 

lL = I-fPLiqL)iqL-0.2lYdqL+ 

+ /' PMiQL-OMfdqL (41) 

J0.54 

As figure (|^) shows this tends to zero in the thermodynamic limit as expected. 



8 Thermalization 

It is very important in a computer simulation to be sure that obtained re- 
sults describe the equilibrium system fluctuations. The main problem in 
the research of this condition is that generally we are interested to the low 
temperature behaviour of the system where the phase space is very com- 
plex. Theorically we should make simulations for a infinite time to be sure 
that equilibrium is joined but this is obviously impraticable. We so make 
everytime errors whit finite-time simulations, but it is possible to prove that 
these ones are comparables with statistical fluctuations. R.N. Bhatt and A. P. 



Young |]11| suggested a thermalization criterion that is so resumable: 



For every disorder realization one define two overlaps 



Qii) = ^T.Siito)S^ito + t) (42) 



and 



Q'it) = ^ E ^^^(^0 + t)sP{to + t) . (43) 

i 

Where to is the esteemed equilibrium time. In the t ^ oo limit the two 
overlaps have the same distribution, but also in a finite-time observation 
(inside an error margin) they are convergent. 
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For small times Q{t) assumes values very near to 1 and the P(q) function 
will have a peak just there. In the same situation Q'{t) will have (starting 
from a random spin distribution) a Gaussian shape peaked around zero and 
with amplitude A^"2. When we increase t, Q{t) and Q'{t) tend to the same 
distribution. We can say to be in equilibrium when Monte carlo statistical 
errors are bigger than finite-time simulation errors. In other words when the 
two overlap distributions are the same function unless statistical errors. We 
report in figure (H) results of this method applied to the lattice with size L=3 
at the temperature T=1.4 (everytime we use tq = to)- 

Considered the originality of dynamic applied in these simulations (we 
have three replicas of the system and so we can't use the Bhatt- Young 
method) , a criterion of thermalization has not yet been defined. However, 
we think that a further confirmation ( thermalizing the free P(q) we have 
equilibrium in all phase space and in particular in the small part we are 
studying with our dynamic) of the joined equilibrium for our system comes 
from the study of the constrained P(q) for different and successive MCS. 
Specifically, the test we made consists in the control, for a given annealing 
schedule, of the behaviour of the distribution function with respectively the 
first third, the second third and the last third of MCS used for the calculation 
of statistical averages. The result relative to size L=5 with 100 samples is 
shown in figure (^. It can be observed that the three curves are pratically 
identical. In fact oscillations between the functions are well inside the er- 
rors from which are afflicted. This test, made on all lattice size, has given, 
with our annealing schedule, good results comforting ourself about the joined 
equilibrium of the system. This is not obviously the faster method ( we are 
trying to find something better) but it is sure and for now this is enough. 

9 The tolerance parameter e 

We have used a particular attention to the behaviour of the constrained P(q) 
varying the tolerance parameter. We have concentred our study in the lattice 
of size L=3 and L=5 performing a lot of simulations in order to understand 
the develop of the overlaps distribution function under the imposition of more 
and more restrictive bounds. 

The work on the lattice of size L=3 is effectued with a small statistic but 
a large range of tolerance parameter values, vice versa for the lattice of size 



12 



L=5 we have used a robust statistic but a little number of e. The two results 
obtained are just compatibles and show, see figures ([T0| )([TT1) , a continuous 
shifting of the tail relative to the negative overlaps towards values more and 
more near to that imposed from the ultrametric bound. This observation 
suggest us to choice the smallest possible tolerance parameter. There is 
obviously some restrictions. Too small tolerance parameter values would 
impose too restrictive bounds in the simulation of the lattice L=3 while too 
big values could be not enough to force the distribution in a decisive way. 
Moreover when the tolerance parameter is bigger than 1 + 2/5qEA we refind 
the free overlaps distribution (see figure (^) and all dynamical spin-fiips are 
accepted from the bound algorhythm. When e decrease the number of spin- 
fiips accepted begin smaller and to have enough statistic we must increase 
the number of MCS for every sample. So we also have to choose e so that 
calculator simulations are not too long in time. After long tests we have 
established that best compromise in our case is e = 0.04. 



10 The constrained P(q) 2" part 

The evidence of first results lead us to deepen the feature of the ultrametric 
structure in the examined system. So we further modify the dynamic without 
changing the philosophy of the approach. Now we fix the overlap ^23 and gis 
to different values (in the first part of the work we fixed gi3 = ^23 — ^Qea)- 
In particular we made our simulations with 

4 1 
^23 - -^Qea qi3 - -^Qea (44) 

or better ^23 = 0.43 and qis = 0.10. By this modification we expect, differ- 
ently from the previous case, that increasing the size of the lattice the third 
overlap qi2 (following the same notation defined in section 5, gi2 = Q) as- 
sumes exclusively the value q = gis how imposed from the ultrametric bound 
so that the constrained distribution function of the overlaps tends to a Dirac 
Delta function peaked around gis. 

The results of these simulations (see also [|l[) are shown in figure (p!2D . It is 
evident in this development of the overlaps distribution function very 

singular. There is a systematic shift, starting from small volumes, either of 
the peak of bigger tallness or the smaller peak toward, probably, the overlap 
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of value gi3 = 0.10. One can see moreover that increasing the lattice size 
the curve tends to narrow and the smaller peak is "absorbed" in the tail of 
the predominant one so that the distribution function assumes the expected 
shape. 

The study of the biggest overlaps develop with the lattice size (see figure 
(p!3D) consents to establish a limit to the shift of the predominant peak of 
constrained P(q). Our best fits with the function 

(45) 

affirm that q!^ax — 0.100 ± 0.028. So in the infinite size limit our function 
will tend to settle about a value compatible with ^13. 

To understand at last which is the shape that the constrained P(q) as- 
sumes in the thermodynamic limit we have effected a study of variation of 
the integral 

h = j]^ PMiQL - Q.l^fdQL (46) 

with the lattice size. Best fits (see figure (0)), always executed by the 
function 

Il = I'- + ^. (47) 

indicate that 

J~ = lim h = 0.000 ± 0.002 , (48) 

SO in the thermodynamic limit the integral (|i6| ) will be zero. In other words, 
being 

^L(gL)(gL - 0.10)2 >o (49) 

also 

Il>0 VL , (50) 



and in order to satisfy the (^Hj) must be verified that 

P{q) = lim Pi(gL) = SiQL - 0.10) (51) 
as required from the ultrametric hypothesis about the states structure. 
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11 Another way to see ultrametric structure 



Another approach we tried to verify the existence of a possible ultrametric 
structure between the low energy states of a realistic spin glasses is suggested 
from the work about the hipercubic lattice . 

We take three replicas Si, S2, S3 of the system with the same disorder 
realization and build the overlaps qi2, qis, 523 • 

In the previous works we fixed by constrained dynamic the value of two 
of three overlaps (at the same Q-^^^ or at differents Q-^^^) and observe the 
behaviour of the third one. 

Now we define an overlaps bound, say [qfi^i, qfix2] (with qfi^i > qfix2), 
and then check (with a not constrained dynamic) whether the largest overlap 
falls inside the established bound. In this case we calculate the difference 
between the other ones 6q = qmid — qmin and plot the distribution function 
P{Sq). We aspect to obtain in the termodynamic limit a Dirac delta function 
peaked around zero. 

We have chosen for our simulations the interval [qEA,'^/5qEA] for lattice 
sizes L=4,6,8 with 400, 200 and 100 samples respectively at the same tem- 
perature T=1.4. 

As we can see in figure (^) increasing the volume of the lattice the curve 
tends to narrow more and more showing the presence of a not trivial structure 
between the low energy states of the model. As well as in the work about 
hypercubic celles, the behaviour of the curve seems to manifest a ultrametric 
feature of the system. 

Our last simulation with all others developed in this work we lead to 
suppose that the spontaneous replica symmetric breaking is a good analitic 
approach for the study of realistics Eduards Anderson Spin Glasses. 



Conclusions 

In this work we describe in a very detailed way the constrained dynamic used 
in Monte carlo simulation of a four dimensional Eduards Anderson Ising Spin 
Glass to understand the nature of the glass phase for these systems. We 
have concentrated our efforts in the study of the low energy states structure. 
After having defined a metric (section 3) in the overlaps space, we have 
effected a long series of simulations (section 7, 8, 9, 10) in 4D lattice of sizes 
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(L=3,4,5,6,7,8) looking for some confirmations about a possible analogy with 
the model SK. From the obtained results, it seems to be possible to affirm 
that a not trivial low energy states structure in the glass phase of a finite 
dimensional Spin Glass really exists. In particular the results shown in figures 
@, ([T2|) and @ claim the concreteness of the hypothesis proposed by G. 



Parisi about the origin and the kind of this structure. The origin should 
be researched in the spontaneous replica symmetry breaking [@] 101^0 and 
the kind of the structure really seems to be the ultrametric one as already 
established for the SK model. 
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Figure 1: Probability distribution of the difference between the two smaller 
overlaps {6q = qmid — Qmin) when the bigger overlap is smaller than 0.4 from 



in a four dimensional Ising Spin Glass for different lattice sizes ||T2 
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Figure 2: Probability distribution of the difference between the two smaller 
overlaps {6q = qmid — Imm) for a fixed larger overlap value q=0.5. The 
temperature is Q.QTcrit.- In the N ^ oo limit, the distribution is expected to 
become a S function at the origin |Tl[ 
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Figure 3: Distribution probability P(gi2) (in y-axes) of starting values that 

g(i2) (in x-axcs) assumes varying the initial spin configurations with 10000 
different starts for sides L=3, 4, 5, 6, 7, 8 in reading order 
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Figure 4: Distribution probability of the overlaps P(q) at T—IA for different 
sizes of the lattice L=3, 4, 5, 6 
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Figure 5: Evolution of the probability distribution of overlaps P(q) with a 
constrained dynamic a T—IA for lattice size L=3, 4, 5, 6, 7, 8. 
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Q max=Q"ii,ax+«*L 
Q"max=0-31±0.09 

c< = 0.85±0.03 
7' = -0.6±0.2 
xV3 = 0. 04590193 
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Figure 6: Fit of biggest overlaps for the constrained P(q) with the lattice 
size by the function J^^^ = /^^^ + in log-log scale 
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I^=I°°+a*L^ 

I"=-0.0001±0.0005 
c< = 0.76±0.03 
y = -2.2±0.04 
xV3 = 0. 4372690 
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Figure 7: Evolution of the integral h = I-f PL{qL){qL - 0.2lf dqi + 
Io.54Pl{ql){ql — 0.54)^dg(L) with the lattice size in log-log scale 
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Figure 8: Evolution of free overlap distributions P(q) and P(q') (in y- 
axes) (see section ^ as a function of thermalization steps (in x-axes). 
to = 10, 10^, 10^, 10^ 2 ■ 10^ 3 ■ 10^ 5 ■ 10^ in reading order 
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Figure 9: The tree curves (basically coinciding in the plot) are for constrained 
P(q) in the first, the second and the third of the run after the annealing 
scheme described in the text), L=5, 100 samples 
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Figure 10: Evolution of the constrained P(q) (in y-axes) as a function of 
overlap (in x-axcs) for different values of tolerance parameter e =0.8, 0.6, 
0.4, 0.3, 0.2, 0.15, 0.12, 0.1, 0.08, 0.06, 0.04. L=3, T=1.4, samples=64. 
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Figure 11: Evolution of constrained P(q) (in y-axes) as a function of overlap 
(in x-axes) for different values of tolerance parameter e — 0.22, 0.12, 0.09, 
0.05, 0.04, 0.02, L=5, T=1.4, samples=100. 
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Figure 12: Study of the behaviour of constrained P(q) with lattice size; L=3, 

4, 5, 6, 7, 8. T=1.4 {qmin Qmax) 
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Figure 13: Fit of the predominant parameter peak of the constrained P(q) 
by the function q^^^ = + ^7 in log-log scale 
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Figure 14: Evolution of the integral II — J^i Pl{(1l){(1l — 0.10)^(igL with the 
lattice size in log-log scale 



29 



P(5q) 



15 




Figure 15: Ultrametric feature of the EA model for lattice of side L=4, 6 ,8. 
We plot the probabihty of the difference between the middle and the smaller 
overlap 5q — Qmid — Qmin when the largest fall inside the range [2/5qEA, Qea] 
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